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ABSTRACT 

Using deep photometric data from WFC@INT and WFI@ES02.2m we measure the 
outer number density profiles of 1 9 stellar clusters located in the inner region of the 
Milky Way halo (within a Galactocentric distance range of 10-30 kpc) in order to assess 
the impact of internal and external dynamical processes on the spatial distribution of 
stars. Adopting power-law fitting templates, with index —7 in the outer region, we 
find that the clusters in our sample can be divided in two groups: a group of massive 
clusters 10 5 Mq) that has relatively flat profiles with 2.5 < 7 < 4 and a group 
of low-mass clusters 10 5 M Q ), with steep profiles (7 > 4) and clear signatures 
of interaction with the Galactic tidal field. We refer to these two groups as 'tidally 
unaffected' and 'tidally affected', respectively. Our results also show a clear trend 
between the slope of the outer parts and the half-mass density of these systems, which 
suggests that the outer density profiles may retain key information on the dominant 
processes driving the dynamical evolution of Globular Clusters. 

Key words: globular clusters: general - methods: observational - stars: Population 
II - techniques: photometric 



1 INTRODUCTION 

Globular clusters (GCs) have been found in nearly 
all galaxies and they have been considered the fos- 
sil records of the formation of their host galaxy (e.g . 
Zinnlll993l; IZepf fc Ashman! Il993l : iBrodie fc Straderl l200rj : 
Marm-Franch et all 120091 ). In t he Milky Way ro ughly 



150 GCs are currently known (|Harrisl Il996l . I2OI0I ) and 
they have a long history of being used in the context 
of stellar evolution (e.g. iBaraffe et alj Il998l ). the struc- 
ture and chemical evolution of the Milky Way (e.g. 
Ilnnanen et alj fl983t ISearle fc Zinn|[l978l . respect ively) and 
the dynamical evolution of collisional systems (|Hut et al.l 
ll992l : lMevlan fc Heggidl 19971 ). 

The radial number density profiles of GCs can be used 
to learn about their evolution becaus e the structu re and evo- 
lution of clusters are closely linked (|Henon| [l961). In recent 
years, much attention has gone to the inner parts of the 
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surface brightness profile to quantify t he number of clus- 
ters that has undergone core collapse (|Diorgovski fc Kind 
1 19861) and to look for the presence of massive black holes 
( Novola fc Gebharddl2006l ). Especially for the densest clus- 
ters, these studies require the high angular resolution of 
space based observatories, such as the Hubble Space Tele- 
scope, to resolve the dense core into individual stars. On 
the other hand, the spatial distribution of the bulk of the 
stars is determined by the interplay between internal and 
external dynamical processes. Internal relaxation processes 
tend to populate the outer regions of the cluster with low- 
mass stars, whereas the external tidal field of the host galaxy 
will strip some of the stars at large radii. Hence, studies of 
the entire profile are useful to gain insight on the overall 
evolution of individual clusters. Since such studies require a 
wide field of view, they are usually done from the ground 
(|Trager et al.lll995l . hereafter TGK95). 

Several self-consistent (static) models for clusters ex- 
ist, usin g different assumptions for the distribution func- 
tion fe.g. |piummedll91ll ; lw"oollev fc Robertsonlll956l ; lMichiel 
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ll963l ; lKing|ll966l ; IWilsonlll975l ). The surface brightness pro- 
file is often used to say something about the correctness of a 
model, but unfortunately it does not always provide enough 
information to constrain the underlying model, i.e. with- 
out the use of additional inform ation such as stellar veloci- 
ties. This was already noted bv lHenonl ljl96ll ) who realised 
that different models could successfully describe the surface 
brightness profile of a GC despite their different underlying 
physical assumptions. 

Most of these models are based on simplifying assump- 
tions, such as spherical symmetry, an isotropic velocity dis- 
tribution, stars of the same mass and a truncation of the 
distribution function. It is, therefore, that empirical tem- 
plates are often used, i.e. simple fitting formulas that do 
not (necessarily) follow from a specific distribution function, 
to fit and compare surface brightness pr ofiles of clusters. 
A well-known example is the iKina (1 19621 ) template that is 
typically used for old GCs and that is characterised by a 
steep drop, or truncation, of the density at the outer edge. 
This radius is called the tidal radius in King's paper, but 
because it does not necessarily coincide with the radius of 
the zero- velocity surface, or Jacobi radius |Baumgardt et all 
l20ld : lKiipper et al.|[201Tjh . we will refer to the radius where 
King model s drop to zero density as the edge radius (r e d ge )- 
lElson et alj (|l987l . hereafter EFF87) used a continuous tem- 
plate (i.e. without a truncation) for relatively young clusters 
in the Large Magellanic Clouds (LMC). The EFF87 tem- 
plate is characterised by a flat core and a power-law enve- 
lope. Th ese are part of the more general family of power-law 
models (|Zhadll99rj ). 

The cluster system of the LMC is a good example of 
how observations of the structure of clusters can be used as 
a proxy for evolution. Because the ages of clusters in the 
LMC span a range from a few Myrs to ~ 12 Gyr, varia- 
tions of structural parameters with age can be studied, and 
when interpre ting this as e voluti on, m uch can be learned. 
For ex ample, lElson et al.l (|l989l ) and iMackev fc Gilmorel 
(|2003bh showed that the core sizes of LMC clusters increase 
(on average) with age, which was also found for clusters 
in the Small Magellanic Clouds (SMC; IMackev fc Gilmorel 
l2003al) and for cluster systems outside the Local Group 
(|Mai'z-Apellanizll200ll ; lBastian et alj|2008r i. Numerical stud- 
ies suggest that this expansion is largely driven by internal 
evolution jMackev et al.l [20081 '). rather than by external ef- 
fects su ch as the details of the orbit around the centre of the 
galaxy (| Wilkinson et al.ll2003h . 

The young SM C/LMC clu sters and young extra- 
Galactic clusters (e.g. Larson 2004) arc well described by the 
continuous power-law models, whereas for old GCs (tidally) 
truncated King models are often used. It seems, therefore, 
logic to conclude that clusters form with extended haloes 
and that a tidal truncation develops over time through in- 
teractions of the cluster with the tidal field of its host 
galaxy. Because the tidal field is strongest at pericentre, 
it is often assumed that this is where the truncation in 
the c luster density is determined {King 19621; Ihmanen et al.l 
1983). Numerical works, however, indicate that the real- 
ity is more complicated. Collision-less iV-body simulations 
of stellar clusters that start with a King-type density dis- 
tribution and that move on eccentric orbits in the Galac- 
tic potential develop an halo of "extra-tida l" stars as they 
shed a fraction of their stellar mass to tides (|johnston et al.l 



Il999l) . Once the unbound material escapes and equilibrium 
is re-established, the outer profile of the cluster does not 
exhibit th e original truncation and instead app roaches a 
power-law (|Oh et al.lll995l ; |Penarrubia et alj|2009h . The sur- 
face density in that halo is shallower than for the bulk 
of the stars and the location where the break occurs is 
an indicatio n of the time elapsed from the last pericen- 
tre passage (|Penarrubia et alTl2009l ). Observational exam- 
ples of this phenomenon have been d i scove red in Pal 5 
Odenkirchen et al.1 120031 ; iDehnen et all |2004). NGC 5466 



Belokurov et al.ll2006h and Pal 14 (|Sollima et ai1l201ll ) 



Tidal stripping is not the only process that can generate 
a break in the surface density profile. Models of star clusters 
that lose stars through two-body relaxation and move on 
a circular orbit, i.e. without the time-dependent tides that 
give rise to tidal stripping, also develop a halo of energet- 
ically unbound stars. These stars were dubbed "potential 
escapers" and they are trapped within the Jacobi surface 
because the criterion for escape is no t only based on en- 
ergy but also on angu lar momentum (jFukushige fc Heggid 
l200d : iBaumgardtlboOll ) . These stars give rise to a very sim- 
ilar break in the outer parts of the surface density profile 
as in the tidal stripping scenario (|Kupper et alj|20ld ). Such 
a break was found i n the surface den sity /brig htness pro- 
files o f the GCs M 92 dTesta et al ll2000t) . Pal 13 dCote et al. 
20021) . Whiting 1 ijCarraro et al.ll2007T ) and AM 4 ijCarraro 
2009h . Often these breaks go together with direct detec- 



tions of the characteristic "S-shape d " tidal features a nd/or 
tidal tails feelokurov et al] l200d : Ijordi fc Grebe! l20ld : 
iNiederste-Ostholt et al.ll2010h . When the details of the orbit 
of the cluster are not available, it is difficult, if not impos- 
sible, to disentangle the relative contributions of tidal strip- 
ping and two-body relaxation to the formation of the break 
in the surface brightness profile. 



iMcLaughlin fc van der Marell (|2005l . hereafter MvM05) 
fit different models to the surface brightness profiles of Milky 
Way GCs using the data of T GK95 and they find that the 
more extended IWilsonl (|l975T ) and EFF87 models provide 
equally good and sometim es even better descriptions than 
the traditional iKind (|l96d ) models. In particular, they show 
that goodness of fit parameters indicate a preference for 
more extended models whenever more data in the outer 
parts are available. They refrain from a detailed interpre- 
tation and conclude that age is not the only parameter that 
determines whether a King-type (i.e. truncated) model pro- 
vides a good description of the surface brightness profile, or 
not. In this study we search for the additional ingredient 
driving the shape of the clusters density profiles. 



This paper is organized as follows. In § 2 we present 
the sample of Galactic GCs included in this study and a de- 
scription of the observations needed for our purposes. In § 3, 
we describe how we obtained the number density profiles on 
which this study is based and the fitting technique adopted 
to derive their structural parameters. In § 4 the number 
density profiles are presented and are interpreted in § 5 in 
connection with possible external and internal processes. In 
5 6 we draw our conclusions. 
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2 THE SAMPLE 

The clusters included in this project are part of a study 
that focu sed on the search for tidal stre ams around Galac- 
tic GCs (|Marti'nez-DeIgado et aljfeoiMl ). For this purpose, 
we have observed a sample of 19 GCs in the Galactic halo 
between 10 and 30 kpc in Galactocentric distance, sampling 
~ 56% of the total number of Galactic GCs in this distance 
range. Since photometry at low Galactic latitudes is severely 
hampered by the presence of field stars, we have excluded 
all those objects located within 20 degrees from the Galactic 
plane with the exception of NGC 2298 and Rup 106, which 
have been observed in spite of their low Galactic latitude 
because of their possible association with the Canis M ayor 
stream (|Bellazzini et all 12004 iForbes fc BridgesH20irj| ). Be- 
cause all the GCs in our sample lie in a 20 kpc wide distance 
range, we can study the effect of a similar external Galactic 
tidal field on the structure of clusters with different proper- 
ties. 

To place the sample in the context of the Milky Way GC 
population, we show their half- mass densities, ph, against 
the Galactocentric distance, Rg, together with the values 
for the rest of the GCs (Figure [1]). A constant mass-to-light 
ratio of 2 was used for this figure and the half-mass den- 
sity is defined as ph = 8M/(37rrf[). Here M is the mass 
and rh is the three-dimensional half-mass radius which is 
estimated from the half-light radius in projection, or ef- 
fective radius (r e g) by c orrecting for the effect of projec- 
tion, i.e. rh = (4/3)r e fi (|Spitzerlll987l ). This is a powerful 
diagram to study the relation between inter nal and exter- 
nal e f fects who both aff ect cluster properties (|lnnanen et al.l 
1 19831 ; I Gieles et al.ll201lh . The lines are 'cluster isochrone re- 
lations' at an age of a Hubble time for clusters with differ- 
ent rnasses_that evolve in a steady tidal field (Appendix B 
in iGieles et af1l201ll ). It assumes that two-body relaxation 
is the dominant process, which is a valid assumption for 
most clusters outside ~ 10 kpc (jGnedin fc Ostriker 1997; 
IVesperini fc Heggiei IhWtT ) . IGieles et all ( |201lf > showed that 



clusters that form deeply embedded within their Roche-lobe 
spend roughly the first half of their lives in an 'expansion- 
dominated' phase and during the second half of their lives 
they are in an 'evaporation-dominated' phase. These authors 
considered relaxation driven evaporation of clusters on cir- 
cular orbits. Because the evolution of clusters on eccentric 
orbits is affected by an additional tidal effect, i.e. tidal strip- 
ping, we introduce a more general reference to these two 
regimes, namely 'tidally unaffected' and 'tid ally affected', 
respectively. According to IGieles et al.l l|201ll ). clusters that 
satisfy 

M<1O 5 M 1^ (1) 
Kg 

are in the tidally affected regime. This constraint is satisfied 
for ~ 25 % of the clusters in our sample: Rup 106, NGC 5053, 
NGC 5466, Pal 5 and NGC 7492. The remaining clusters of 
our sample will be referred to as tidally unaffected clusters. 

To derive the mass of the clusters (Table 2), we have 
taken the V-band mass-to-light ratios found by MvM05 and 
the luminosities derived from the integrated V magnitudes 
of our profile^. Moreover, trying to complete our analysis 
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Figure 1. Half- mass densities against Galactocentric distances 
of the 141 globular clusters (black dots) of the Harris catalogue 
(Harris 1996). The clusters used in this study arc highlighted with 
a (red) star. All symbol sizes scale with the mass of the cluster. 
Cluster isochrone relations from Gieles et al. (2011) are over- 
plotted for clusters with different masses. Clusters that align with 
the horizontal lines are in the tidally unaffected regime, whereas 
clusters that align the diagonal line are in the tidally affected 
regime. The exact separation between these regimes is given by 
equation JlJ. 



of the dynamical state of each cluster, we have used the 
expression 



t T h — — 



0.138 



Mri 



1/2 



Myr 



(2) 



mln(0.11M/m) ^ G 

bv lSpitzer fc Hard l|l97lh to obtain the half-mass relax- 
ation time in Myr. In this formula, m is the mean stellar 
mass that has been set to 0.5 Mm. The constant of 0.11 was 



found bv lGiersz fc Heggiei (|l994T l for clusters of equal mass 
stars. 

The six-dimensional phase space information of a GC 
is essential to derive its orbit in the Galaxy and evaluate its 
interaction with the Milky Way tidal field. While accurate 
positions on the sky, distances and radial velocities are avail- 
able for all the cluster of our sample, proper motions are dif- 
ficult to obtain for these distant objects. We only have found 
orbital param e ters f or 10 of our GCs in the compilation by 
iDinescu et al.1 (|l999T), Absolute ages included in Table 1 are 
extracted from the database by Forbes fc Bridges! l|2010T ) ex- 
cept for NGC 6229 jBorissova et alJll999h .' Coordinates and 
dista nces are taken from the Harris catalogue (|Harrislll996l . 
l2010h . 

2.1 Observations 

As already reported in the previous Section, the observa- 
tional data used in this paper come from a photometric cam- 



1 To convert V magnitudes into luminosities we adopted the dis- 



tance moduli and reddenin g from Harris 1 2 010T) and My q = 4.83 
iBinnev fc Merrifieljfl998h 
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Table 1. Sample of Galactic GCs included in this project. Coordinates and distances were obtained from the updated version of Harris 
catalogue (Harris 1996, 2010) while radii have been inferred from the structural parameters published in MVM05, except for Rup 106, 
taken from Harris. Absolute ages are those computed by Forbes & Bridges (2010) except for NGC 6229 (Borissova et al. 1999) and 
Galactocentric distances during apo- and pericentre of the clusters are obtained from Dinescu et al. (1999). 



paign aimed at searching for tidal streams around Galactic 
GCs. Due to the low surface brightness of known Galactic 
tidal streams (~32 mag/arcmin 2 ), a suitable combination 
of wide field and exposure time was required. We have used 
the Wide Field Camera (WFC) located at the Isaac New- 
ton Telescope at the Roque de los Muchachos Observatory 
(La Palma, Spain) and the Wide Field Imager (WFI) at the 
ES0 2.2m telescope at La Silla Observatory (Chile). The 
WFC covers 34x34arcmin with 4 identical chips while the 
WFI has a similar field of view (FOV) of 34 x 33 arcmin 
distributed in 8 chips. 

For those clusters with angular sizes comparable to the 
FOVs of the instruments, it was necessary to include an 
additional field to cover the outer regions, where we focus 
on. The total exposure times in B and R bands were 4x900 s 
and 6x600 s respectively reaching a limiting magnitude of ~ 
5 mag below the turn-off in the colour-magnitude diagram 
(CMD). With such deep photometry we have been able to 
detect cluster members in the outer parts of these systems 
where the stellar density is critically low. 

The images were processed using the imred package 
and standard routines in IRAF. Photometry was obtained 
with the PSF- fitting algorithm of DAOPHOT II/ALLSTAR 
l|Stetson|[l987h . One of the main advantages of using this 
software is that it provides criteria to reject extended objects 
as background galaxies so we only included stellar-shaped 
objects in our final photometric catalogue. This decision is 
taken based on the sharpness parameter for which we set 
the limitation |sh|^ 0.5. During each observing run a set of 
photo metric standard stars from the catalogue bv lLandoltl 
(|1992|) have been observed. They have been used to derive 
the transformation between the instrumental magnitude and 
the standard Johnson system. We adopted the atmospheric 
extinction coefficients provided by both observatories. 



3 NUMBER DENSITY PROFILES. FITTING 

To construct the radial density profiles, we selected the 
bona-fide cluster members along the main-sequence (MS). 
First, we fitted the CMD o f each cluster with a set of 
isochrones from Marigo et al] (l2008h wit h suitable age and 
metallicity from lForbes &i Bridges! (120101) and adopting the 
distance and reddening bv lliarrisl l)201fj ). We defined a box 
in the CMD containing all the stars with (B-R)i so -(B- 
#)cmd| < 0.15 and B to < B < B m „, where (B - R) iso 
and (B — R) C md are the isochrone and stars colours, respec- 
tively, B to the turn-off magnitude and -B max a limiting mag- 
nitude defined to be 2.23 ^ -Bmax — Bto ^ 4.48 depending 
on the cluster. The adopted magnitude range represents a 
good compromise to maximize the number of cluster objects 
ensuring a good level of completeness (0 ^ 90%) estimated 
from artificial star experiments in the most crowded region 
used in the analysis for a s ubsample of relatively dense GCs 
(see iBellazzini et al.ll2002l ). The comparison with the theo- 
retical isochrones indicates that in all the GCs of our sam- 
ple, the stars that satisfy the adopted selection criterion lie 
in a similar mass range, between 0.77 ^ M max ^ 0.86 M© 
and 0.52 < M min < 0.69 M e , where M max and M min are 
the maximum and minimum stellar masses found in the MS 
of these clusters. Although the adopted selection criteria in 
both sharpness and location in the CMD minimize the con- 
tamination from back/foreground field stars and galaxies, 
the presence of some intruders in the sample is unavoidable. 
The fraction of contaminating objects has been estimated 
for a subsample of clusters for which suitable control fields 
located at the same Galactic latitude but ~ 3 deg away from 
the cluster position. We estimated a fraction of outliers sat- 
isfiying the selection criteria described above of 0.1-0.3%. 
Morever, such objects are expected to be homogeneously 
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distributed across the field of view and should therefore not 
alter significantly the profile shape. So, we can can safely 
consider negligible the effect of such intruders in the derived 
profiles. 

We counted all those stars included in our box that are 
contained in concentric annuli, centred in the cluster cen- 
tre coordinates with fixed width in logarithmic scale. The 
number of counts per unit area has been obtained by divid- 
ing the number of stars that satisfy the above criteria by 
the corresponding area covered by the annuli. The area of 
each annulus has been corrected for the region not covered 
by our observations (gaps between the chips, borders trun- 
cation, ect). The error on the density has been estimated 
using the standard error propagation formula and assuming 
a Poisson statistic for star counts. 

For most of the clusters of our sample the half-mass 
relaxation time is smaller than their age. It is therefore im- 
portant to investigate the effect that mass segregation could 
have on the resulting profiles by splitting the MS in two sub- 
samples of equal magnitude extent covering different mag- 
nitude ranges. We found no signs of significant differences in 
the obtained density profiles for all the clusters of our sample 
along the entire cluster extent, with An(r) ^ 1.5% (see Fig- 
ure for the illustrative case of NGC 1904). The Kolmorov- 
Smirnov test indicates a probability of ~ 98% that the radial 
profiles obtained for the two different subsamples of MS stars 
are extracted from the same population. So we can neglect 
the effects that mass segregation could have on our derived 
profiles. 

The long exposure times needed to reach a lower limit- 
ing magnitude while allowing to efficiently study the outer 
region of these GCs, causes severe incompleteness in the 
inner arcminutes which prevent us from studying the struc- 
ture in the innermost portion of the cluster. As a solution, we 
have used the catalogue of surface brightness profiles of GCs 
in the Milky Way published by TGK95 (with the exception 
of Rup 106, not included in that catalogue). For each clus- 
ter we have converted the TGK95 profiles into star counts 
adopting the relation 



n(r) = C-0.4 nv{r) 



(3) 



Then we fitted both our and the TGK95 dataset simulta- 
neously, leaving the vertical scaling factor between the two 
profiles, C, as a free parameter in the fitting model. 

In the fitting procedure we optimized the combined 
chi-square (x 2 ) of both datasets. When fitting, we noticed 
that certain profiles in the TGK95 dataset have a fraction 
of points with very small (probably unrealistic) error-bars, 
which would completely drive the fit. Furthermore, in sev- 
eral cases, the data from TGK95 in the very inner parts of 
the clusters show some small scale structure in the profiles, 
which can be either intrinsic or due to systematic effects. 
In order for the fit to be determined mainly by the over- 
all profile shape rather than a few points/structures in the 
TGK95 data, we decided to limit the smallest error-bar of 
the TGK95 dataset to 0.1 dex in density. 

The fitting model for the profiles was: 



n(r) = n B G + n c i(r) 



(4) 



where tibg is the density of background/foreground 
sources and n c i(r) is the model for the cluster density pro- 
file. In this paper we considered two different models. As in 
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Figure 2. Example of MS stars selection to calculate the number 
density profile of the cluster NGC 1904. The lower panel shows 
the profiles for different masses ranges (0.68 C M ^ 0.83 Mq, 
open red circles; 0.54 ^ M < 0.68 Mq, solid blue circles). 



EFF87, a power-law model with parameter 7 controlling the 
shape of the profile : 

1 



nci(r) = n 



(l + ir/ro) 2 )^ 2 



(5) 



and the King (1996) profile, parametrized by the Wo 
parameter: 

n c i (r) = n x F(r,r , W ) (6) 

Because I Kind (|l966t ) models do not have an analytical 
solution for the surface density, a grid of models was pre- 
computed by numerical solving the Poisson equation. 

In order to fit the density profiles, instead of standard 
gradient descent methods (e.g. Levenberg-Marquardt), we 
decided to employ a Markov Chain Monte-Carlo (MCMC) 
technique to better explore the possible covariance between 
parameters. For each cluster, we ran multiple MCMC chains, 
and then from all the chains we picked up the values of 
the parameters giving a best fit to the data, while the ID 
posterior distributions of the parameters gave us estimates of 
the error-bars, which were often large due to the covariance 
between them. Since our data combined with the TGK95 
data covers very large range of radii, for a few clusters we 
do see systematic deviations of the observed profiles from 
our simple models. In these cases obviously the error-bars 
on the parameters are only indicative. In a few cases the 
background is not well sampled by our density profiles. In 
our fitting procedure, the sky level is a free parameter with a 
fiat prior between and the average of a few outer points of 
the density profile. This parameter is sampled by the MCMC 
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chain and the values of all the free parameters have been 
determined from marginalized ID posteriors, so the errors 
on the structural parameters are in some cases significantly 
influenced by the errors of the sky determination. 



4 RESULTS 

The derived number density profiles are shown in Figures [3] 
to [T] In order to obtain a better vision of the divergences be- 
tween the King (blue dashed line) and power-law (red solid 
line) fittings, we have included the profiles both in logarith- 
mic and linear scales in the left and right panel, respectively. 
We have also indicated the start point of our data with an 
arrow pointing down, whereas upward pointing arrow indi- 
cates the last data point of the TGK95 profile for each clus- 
ter (with the exception of Rup 106 which is not contained in 
that study). 

The tidal radii previously estimated by MvM05 for 
these clusters are marked in the profiles as vertical dashed 
lines. For most of the cases, the known tidal radius does 
not seem to be compatible with the profiles obtained in this 
work. In clear cases as NGC 1851, NGC5824 or NGC 6229, 
the distance at which the truncation on the profile is ex- 
pected using the TGK95 profiles alone seems to be located 
within the continuous power-law profile found in our data 
set. This fact shows the risk of computing these parameters 
from FOV-limited photometry. 



Since in our power-law model n c \(r) 



for large 



distances from the cluster centre, the parameter 7 represents 
the slope of the outer regions of the profiles and we can use 
it as an indicator of the overall shape of the profile. To have 
a better vision of the variation of the number density profiles 
in connection with the physical properties contained in Table 
2, they have been sorted in ascending order of 7. 

It is visible both in the density profiles and the infor- 
mation shown in Table 2, that the tidally affected clusters 
NGC 7492, NGC 5466, Pal 5, NGC 5053 and Rup 106 are the 
low-density subsample with steeper profiles. The effects that 
variable tides have in the overall structure of these clus- 
ters is evident with the appearance of breaks in the outer 
parts profiles in some of the clusters, with NGC 5466 and 
Pal 5 as best examples. We also know of the existence of 
more or less coherent tidal tails e merging from NGC 505 3 
(lLauchneret al.ll2006ri, NGC 5 466 |Belokurov et al.ll2006h. 
Pal 5 (|Odenkirchen et al||200ll ) and NGC 7492 (|Leon etall 
I2000T I. In this context, the shallower profile observed in the 
outer regions of NGC 5466, Pal 5 and possibly NGC 7492 is 
generated by the dramatic impact that tidal stripping has in 
populating the outskirts of the cluster. Given that the sec- 
ond slope beyond the break only modifies slightly the global 
profile, we proceed with the analysis of these clusters in the 
same conditions as the rest of the sample. 

On the other hand, massive clusters and more generally 
denser systems, present a relatively flatter profile with a 
remarkable continuous power-law distribution. This group 
coincides with what we here refer to as tidally unaffected 
clusters and in these cases we do not detect any signal of 
tidal truncation on the profiles or the existence of potential 
tidal debris as in the low-density group. 

Using the 7 parameter as an indicator of the overall 
structure, we conclude that, as a general rule, tidally affected 
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Figure 8. King model fit for NGC 1851. The fit made using only 
TGK95 data (black dots) and the combined sample constructed 
adding the densities measured in this work (open dots) arc shown 
with red solid and blue dashed lines respectively. Vertical line 
indicates the position of the edge radius derived by MvM05. 



clusters present 7 > 4 while the tidally unaffected group is 
confined to the range 2.5 < 7 < 4. This a grees with the 
profile obtained for an isolated M31 GC by Mackev et al.l 
|201(]j), which ranges from 7= -2.5 in the inner parts to 
7= -3.5 far from the cluster centre. It is also interest- 
ing to note the coincid ence with the result presented by 
iBaumgardt et all (|2010T l. where the existence of two pop- 
ulations in the Galactic GCs was found, based on the ra- 
tio of their half-mass to Jacobi radii. Our tidally affected 
(tidally unaffected) clusters are referred to as tidally filling 
(compact) in that work. 

We want to remark the problems that King model has 
to describe the profiles of NGC 1851, NGC 5272, NGC 5824 
and NGC 7078, four of the most massive clusters included in 
our sample, as well as NGC 5466, NGC 5694 or NGC 6229. In 
the case of NGC 7078 it is not surprising since (as NGC 1904) 
is clas sified as core-collapsed GC in the iDiorgovski fc Kind 
(|l986h catalogue. In all the other cases, the inclusion of our 
data in the outermost part of the clusters drives the fit to- 
wards small core radii and the failure of the King model in 
describing the overall cluster structure (see Figure [8} . On 
the other hand, we find NGC 5634 that is clearly difficult to 
fit with a single power-law. 

Structural parameters derived from the fitting as well 
the masses and half-mass densities inferred from our pro- 
files (using a power-law description) have been included in 
Table 2. The uncertainties both in Wo and 7 show how im- 
portant is to have the inner regions to constraint the overall 
results: in the extreme case, Rup 106, the lack of informa- 
tion of the central parts of the cluster generates problems 
for both models to obtain a reasonable fit. 



5 DISCUSSION 

5.1 Power-law template vs. King model 

The number density profiles presented in this paper show 
that, in most of the cases, the outer parts of the clusters ex- 
tend to large distances from the centre, generally beyond 
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Figure 3. Number density profiles for NGC1261, NGC 1851, NGC 1904 and NGC2298. Blue dashed line corresponds with the best 
King model fitting while red line shows the best power-law fitting following the expression described in Section 3. Upper arrow shows the 
initial data point we have obtained from our data while lower arrow represents the last TGK95 data point. Vertical dashed line indicates 
the tidal radius obtained by MvM05. Note that the vertical range of the plots might differ from cluster to cluster. 



the edge radius inferred from the previous studies based 
on datasets restricted to the innermost region. As a con- 
sequence, the systematic application of the King model to 
previous data and the lack of complete radial profiles seems 
to have generated the idea that these profiles are truncated, 
whereas here we show that in some cases the profile is well 
described by a continuous power-law down to our detection 
limit. 



To analyse the effect that increasing the observed FOV 
has on the derived properties of GCs, we have plotted our 



best fit 7 and r e d gc against those derived by MvMOE0 (which 
are based on the TGK95 surface brightness profiles) in Fig- 
ure [9] (upper panels). From this comparison, the presence of 
clear trends is evident: for tidally unaffected GCs (7^4), we 
find slightly steeper profiles (larger 7), while for the tidally 
affected clusters NGC 5466, Pal 5 and NGC 7492 the profiles 



2 Note that MvM05 use 7 to describe the power-law decline of the 
3- dimensional number density profile. We, therefore, subtracted 
1 from their values to be able to compare to our 7, that describes 
the decline of the number density profile in projection, i.e. 7313 = 
72D + 1. 
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Figure 4. Same as Fig.|3]for NGC4147, Rup 106, NGC4590 and NGC 5024. 



are slightly shallower (smaller 7) . This result is expected be- 
cause of the presence of tidal debris in the outermost parts 
of these clusters which were not contained in the original 
TGK95 data. New edge radii are also remarkably bigger 
than the previous values reported by MvM05 with a mean 
variation of ~+40 %, excluding the 4 GCs which present 
smaller derived edge radii and those which are poorly fit- 
ted by King models (see below). With these new results, we 
do not detect signatures, at least in the tidally unaffected 
group, of truncation on the profiles. 

To investigate the suitability of the King model and a 
power-law to describe our GCs, we compared the computed 
X 2 ■ In Figure [5] we plot the rate of the \ 2 °f t ne power-law 
fit, Xpli over the \ 2 following from the King model fit, Xk> 
as a function of 7 (bottom-left panel) and the same ratio if 



we only include the TGK95 data (bottom-right panel). Ac- 
cording to these results, the power-law template is a better 
approximation to the observed number density profiles for 
~ 2/3 of the clusters in our sample with the remarkable ex- 
ceptions of NGC 5634 with Xpl/Xk ^ 3. Again, the density 
profiles of the tidally affected clusters (with the exception 
of NGC 5053) are better reproduced by power-law fits. The 
same result can be found also using only the inner part of 
the profile from TGK95. It is interesting to note the striking 
difference in the Xpl/Xk ratio of NGC 6864 (member of the 
tidally unaffected group of clusters) which is better fitted 
by a King model when only the inner region is considered, 
while showing a power-law shape when the entire profile is 
considered. 

The reason for the deviation from the King model pro- 
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Figure 5. Same as Fig. El for NGC 5053, NGC5272, NGC 5466 and NGC5634. 



file in the group of tidally affe cted GCs c an be found in 
the working hypothesis made in iKind l|l966f ) models on the 
radial truncation at the zero-energy surface. In fact, while 
stars with energy above the critical energy (i.e. a velocity 
larger than the local escape velocity) are actually expected 
to evaporate from the cluster, in real clusters they remain 
marginally bound to the system for a timescale compara- 
ble to the cluster orb ital period arou nd the parent galaxy 
(I Lee fc Ostrikeilll987r ). or even longer ()Fukushige fc Heggid 
l200d ). This effect produces a significant overdensity of stars 
in the outermost regions of the cluster with respect to the 
predictio n of a King model with out affecting the inner profile 
(see also I Johnston et al.lll999T ). This effect is more evident 
in clusters in the evaporation regime (i.e. tidally affected) 
and those subject to a stronger interaction with the exter- 



nal tidal field. The assumption of the King edge radius as a 
real physical limit of the cluster, leads to the use of the term 
'extra-tidal stars' for those cluster members outside this ra- 
dius. This has been confir med in several GCs d uring the 
last years, with NG C 1851 dOlszewski et al.ll2009|V IC4 499 
I Walker et al.ll201lh or NGC 5694 jCorrenti et al]|201ll ) as 
recent examples, and can be noticed also in an important 
part of our profiles (see Figures [3] to [7} . Of course, despite 
of the above limitations, King models remain a valid repre- 
sentation of a GC (at least in the inner part of the profile). 
However, because the power-law template is a reasonable 
choice to describe the derived number density profiles, we 
decided to use the outer slope 7 as an indicator of the over- 
all structure of our target GCs in the following sections. The 
main consequence of assuming a power-law as the best tern- 
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Figure 6. Same as Fig.|3]for NGC 5694, NGC5824, Pal 5 and NGC 6229. 



plate to model a radial profile is that GCs have an infinite 
extension. Of course we assume that at some point there ex- 
ists a physical limit at the first Lagrangian point or Jacobi 
radius. 

In the following we analyse the impact that both ex- 
ternal and internal factors have in the observed values of 
7- 



5.2 External factors 

It is well known that the structure of a GC is in- 
fluenced by the interaction with the host Galactic 
potential (Ostrikcr Retail Il972l ; iGnedin fc Ostrikerl 1 19991 



Ijohnston et al.l 1999h . In fact, at every passage across the 
Galactic disc and at every pericentric passage, a transfer of 



kinetic energy occurs in the form of compressive (disc) and 
tidal (bulge) shocks. As a result, the escape of high veloc- 
ity stars is triggered while the overall energetic budget of 
the cluster is continuously altered. The theory of the struc- 
tural evolution of GCs as a result of external tidal effects 
has been extensively investigated with analytical approxi- 
mations JChernoff fc Shapirdll987|). t hrough Fokker-Planck 
modelli ng dGnedin fc Ostrikei 1 19991 ) and TV-body simu- 
lation s (|Baumgardt fc Making 120031 ; iGieles fe Baumgardtl 
|200S| ). 

It is therefore interesting to check how the strength of 
these external factors correlates with the shape parameters 
of our sample of clusters. This task requires the knowledge of 
the cluster orbits to quantify the impact of the Galactic tidal 
field. Unfortunately, among the 19 clusters of our sample, 



Outer density profiles of Galactic globular clusters 11 




Figure 7. Same as Fig. [3 for NGC6864, NGC7078 and NGC 7492 . 



only 10 of them have known proper motions so that orbital 
parameters can be derived only for this subsample. 

As a first test we estimated the relative importance 
of the destruction rate due to relaxation and shocks. For 
this purpose we adopted the d estruction rates du e to disc 
+ bulge shocking computed by lAUen et al.l i|200d ) and the 
contribu tion of relaxation take n from the iV-body simula- 
tions of iBaumgardt fc Makind (|200:# l. In Figure \W\ the 
ratio t'disc+buige/^rix is shown for the 10 GCs of our sam- 
ples with known orbits. For two clusters shocks turn out to 
be important (Pal 5 and NGC 5466) , both belonging to the 
tidally affected group. For the other GCs (all belonging to 
the tidally unaffected group) evaporation due to relaxation 
is the dominant effect. The position of these clusters in the 
7 — fdisc+tmige plane is shown in the left panel of Figure 1111 
It is evident that the two tidally affected clusters, character- 
ized by steeper profiles, are subject to significantly stronger 
shocks with respect to the group of tidally unaffected clus- 
ters. 



3 As noticed bv lGnedin fc Ostrikerl jl997f> , the combined effect of 
shocks and relaxation is generally non-linear. Therefore, although 
a rigorous decoupling of these two effects is not possible, such an 
approach must be considered an approximated way to quantify 
the relative importance of the two effects. 
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Figure 10. Ratio between the disc+bulge shocking (from Allen 
et al. 2006) and relaxation from (Baumgardt & Makino 2003) 
destruction rates as a function of 7. The two clusters belonging 
to the tidally affected group (Pal 5 and NGC 5466) are indicated. 



In absence of accurate orbital parameters it is not pos- 
sible to check the above correlations for the entire sample 
of GCs. However, as a general rule, it is expected that the 
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Figure 11. Left: 7 parameter as a function of the disc+bulgc shocking destruction rate v from Allen et al. (2006). The two clusters 
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Table 2. Physical and structural parameters. Masses, densities and half-mass radius relaxation times derived as described in Section 2 
and using the new estimations obtained from our best power-law fitting. Ranges of confidence for Wo and 7 have been included to 
show the discrepeances in some of the clusters, specially Rup 160. 



interaction between low-density clusters and the Galaxy be- 
come more important because their densities are close to 
the background density of the Galaxy. In the right panel of 
Figure [TT] we compare the half-mass radius density with the 
density of field stars within the Jacobi radius calculated as 

PJ ~ 5.376 ^y 2 M eP c- 3 (7) 

according to appendix B in iGieles et all (|201ll ). As ex- 
pected, the tidally affected clusters are restricted to the 
log(ph/pj) < 3 region of that plot while the rest of 7 < 4 
clusters are at least 3 orders of magnitude denser than p,j. If 
we add to our sample those included in MvM05 (grey circles 
in Figure [TTT) we find that clusters in the range log(ph/pj) ^ 
3 follow the same tendency to be represented by 7 < 4. On 
the other hand, there is a group of objects at log(ph/p,j) ~ 1 
and close to 7 = 2 which appear to stray from this behaviour. 
We will discuss the position of this group of clusters in the 
next section. 

As a further test we also correlated the outer 
slope parameter 7 with the main cluster orbital pa- 
rameters total energy, z-component of the angular mo- 
mentum, pericentre radius, orbital period and eccen- 
tricity (E, Lz , R Pe ri , P and e respectively) taken from 
iDinescu et ail (| 19991 ) and constructed an additional param- 
eter F=(Rg— R P eri)/(Rapo— Rpcri) which is an indicator of 
the current position along the orbit of the clusters. With 
this definition, F=l when the cluster is at the apocentre 
while F=0 when it is at its pericentre. In these planes the 
orbits of the two tidally affected clusters are those with the 
largest eccentricity and the smallest modulus of angular mo- 



mentum. This confirms again that eccentric orbits, which 
make the cluster subject to a significant variation of the sur- 
rounding external field, can lead clusters toward the evapo- 
ration regime. A group of tidally unaffected clusters exibit 
orbits with similar eccentricities and angular momentum. It 
is interesting to note that all these clusters have been sus- 
pected to be ac creted (e.g. NGC1851, N GC1904, NGC2298 
and NGC4147; iForbes fc Bridge j|2010h or are very massive 
(NGC7078). Accreted clusters might have different evolu- 
tionary histories from Galactic GCs, passing a significant 
part of their evolution in smaller systems as dwarf galaxies 
before the accretion, so this could explain why accreted and 
in-situ clusters present a different overall structure, even if 
they have similar orbital parameters. Mass represents also 
an important parameter since massive clusters could resist 
the strong tidal interaction that lead towards high gamma 
values. As we have commented, clusters on highly-eccentric 
orbits as NGC 5466 and Pal 5 are influenced by the interac- 
tion with the Milky Way significantly, leading to the forma- 
tion of an external (and different) power-law distribution. 
No other significant trends are visible in the plots shown in 
Figure [T2l between 7 and the other orbital parameters. 



5.3 Internal evolution 

After evaluating the influence of external factors on the ob- 
served number density profiles, we looked at the correlations 
between the general physical parameters of the cluster with 
the slope 7. As in stellar physics, the global mass is always 
one of the first suspects to drive the behaviour of the struc- 
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Figure 12. 7 as a function of the orbital parameters: eccentricity (upper left panel), period (upper right panel), pericentric distance 
(middle left panel), energy (middle right panel), angular momentum (bottom left panel) and phase (bottom right panel). The location 
of the two clusters of the tidally affected group (Pal 5 and NGC 5466) are indicated in each panel with open circles. 
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Figure 13. 7 as a function of mass and density. The dashed line in both panels indicates the tentative border between the two categories 
of clusters (at 7 = 4). Grey points mark the sample of Galactic GCs studied by MvM05. 
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ture of a GC and we have tested this correlation in the left 
panel of Figure [13] 

It has been already noticed (Section 4) that all the clus- 
ters of the tidally unaffected subsample, formed by massive 
clusters, present flatter radial profiles while the less massive 
tidally affected clust ers have large value s of 7. A similar re- 
lation was found by iMcLaughlinl l|200Cl ) using in that case 
the concentration log(r d go /r c ) as a function of the cluster 
luminosity. 

A possible interpretation of this behaviour can be 
searched by looking at the role of binaries in the long- 
term evolution of GCs. Low-mass clusters are able to form 
and retain more binaries than massive clusters where they 
are more efficiently destr oyed in single-binary and binary- 
binary close encounters (|Giersz fc Spurzeml f 19941 : ISollimal 
120081 : iFregeau et all 120091 '). While the natural evolution o f 
GCs is toward high-concentration structures (|Henonl [l 96 lh 
a large fraction of binaries prevent the contraction of the 
core and maintain the cluster in a quasi-steady state of bi- 
nary burning where the continuous loss of energy due to 
evaporation is balanced by the e n ergy heating of binarie s 
(|Gao et alj Il99ll : IVesperinil 1 19981 : IFregeau fe Rasid l2007t ). 
Unfortunately, the binary fraction has only been estimated 
for NGC 4590 (tidally u naffected) , NGC 505 3 and NGC 5466 
(tidally affected) (see ISollima et alj [20071 ) and appear to 
cover the same range (9.5 < £ < 14.2%) within the errors. 
Although the available data do not allow to confirm the sce- 
nario proposed above, other studies have suggested an anti- 
correlation between the mass and the fra ction of binary sys- 
tems in larger samples of Galact ic GCs (|Milone et al .1 [20081 ; 
ISollim a 2008: Sollima ct al. 2010|) ■ I n addition, mass loss due 
to stellar evolution drives the dynamical evolution in early 
stages and causes low-mass stellar systems to expand to their 
tidal boundary fas ter than their more massive counter-parts 
jGieles et al.ll20l"oT ) . This, together with the external factors 
discussed in the previous section, could explain the tendency 
found in the low-mass tidally affected clusters to have promi- 
nent cores, large half-mass radii and larger values of 7. 

We have to analyse the exception of NGC 2298 and 
NGC 4147 in the 7 vs. log(M) plot. Despite their low 
masses, these clusters are included in the tidally unaf- 
fected group according to Table 2. Special birth conditions 
or an external origin could explain these differences. This 
last hypothesis has been widely supported fo r NGC 4147 
([Bellazzini et al.ll2003l ; iForbes fe Bridges! 120101 ) since it lies 
in the projected position path of the Sagittarius tidal stream 
sharing the same energy and angular momentum. Similarly, 
NGC 2298 has been associated to the controversial Canis 
Majo r dwarf galaxy (|Martin et al.l 120041 : IForbes fe Bridges! 
|2010| ) although it presents a retrograde orbit incompatible 
with that of this over-density but also with the general rota- 
tional pattern of the Milky Way. In this picture, GCs formed 
in less denser systems, might live as isolated (tidally unaf- 
fected) clusters before the accretion by the Milky Way. 

The derived outer slope 7 vs. the half-mass density ph 
is shown on the right panel of Figure [131 Also in this case, a 
striking correlation is evident in the sense that denser clus- 
ters have on average flatter profiles. In this case, denser clus- 
ters are both expected to be more resistant to the external 
tidal stress exerted by the Milky Way (see Section 5.2) and 
dominated by the effect of relaxation (eq.[2}. In Figure [T3l we 
have also included the sample of clusters studied by MvM05. 
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Figure 14. 7 as a function of age/t,.^. Grey dots mark the Galac- 
tic GCs of the MvdM sample. 



The trend found in the relations between mass, half mass 
density and 7 are also visible here although the presence 
of a branch of massive clusters with 7 ~ 2 and reaching 
very low densities decreases the significance of the 7 — ph 
correlation. 

Of course, the MvM05 sample contains an order of mag- 
nitude larger number of objects and constitutes a more ro- 
bust statistical sample to check such correlations. However, 
it is noticeable that for most of these object the TGK95 
profiles cover a very limited FOV (including only the first 
central arcminutes). As shown in Section 4, an incomplete 
radial coverage can significantly alter the estimate of the 
profile slope leading to a more uncertain determination of 
both 7 and ph- It is therefore possible that the surprisingly 
low half-mass density of these object is due to an underes- 
timate of 7 and, consequently, ph- 

Connecting the structure of our sample of Galactic GCs 
with their dynamical age is not simple because of t heir nar- 
row range of relative ages (|Marin-Franch et al . 2009]). In Fig- 
ure [2J we show the parameter 7 as a function of the ratio 
age/t r h- On average, clusters with lower ratios of age/trh 
have larger values of 7 even if the large scatter prevents 
from any firm conclusion on the existence of any reliable 
correlation. The same behaviour is visible among the clus- 
ters in the MvdM05 sample with the exception of the group 
of globulars with 7 ~ 2 already mentioned above which show 
a lower age/t T h- It is difficult to say something about the re- 
lation between 7 and the dynamical age of clusters. This 
is firstly because clusters in the tidally unaffected regime 
have expanded such that their relaxation times have become 
a fixed fraction of their age, roughly equal to ~ 1/1 for 
clusters of equal masses ([Henonlll965l : lGoodman|[l98i ) and 
~ 1/3 for clusters w ith a globular type stellar mass- function 
|Gieles et alj|201ll ). The majority of clusters in our sample 
is in this regime and forms a cloud of points in Fig. 1141 Sec- 
ondly, the evolution of clusters with an age roughly equal to 
t r h is probably driven by external factors due to their low 
density, hence t T h is not telling us much about the evolution 
for these objects. 
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6 CONCLUSIONS 

We have presented the number density profiles of 19 Galactic 
GCs obtained with deep wide-field photometry. The compar- 
ison of our results, fitting both King models and power-law 
templates, with those derived by previous studies shows the 
importance that data with a wider FOVs have on our under- 
standing of these stellar systems. In particular, the edge radii 
estimated from this analysis have been found to be ~ 20% 
larger than what found in previous works based on shal- 
lower surface brightness profile. The number density pro- 
files of our sample of GCs can be reasonably fitted by both 
King models and power-law templates. The latter ones ap- 
pear to be a better representation for ~ 2/3 of the observed 
profiles including some of those cl usters showing evidences 
of tidal tails (Pal 5 an d NGC 5466, lOdenkirchen et al.ll200ll ; 
iBelokurov et alj |2006|). in agreement with the prediction of 
JV-body simulationsjLee et aljfeOQfj ; iPenarrubia et al.ll2009l ; 
I Johnston et alJll999f T~ 

It has been found that the slope 7 is linked to the mem- 
bership of the cluster to the group of " tidally affected" or 
"tidally unaffected" clusters (defined bv lGieles et alj (|201lh 
on the basis of the position in the ph(M, Rg) plane). The 
group of tidally affected clusters, constituted mainly of 
low-mass objects with prominent cores, is characterized by 
steeper profiles (7 > 4) while tidally unaffected clusters show 
flatter profiles extending to large distances from their com- 
pact cores. 

We investigated the dependence of the outer slope of 
the density profile 7 on the internal structural evolution of 
the cluster (based on the relaxation process) and on ex- 
ternal factors (i.e. tidal shocks) using a subsample of 10 
GCs with known orbits. External factors have been found 
to be dominant in the two tidally affected clusters (Pal 5 and 
NGC 5466) , which have preferentially more eccentric orbits 
and larger destruction rates than tidally unaffected clusters. 
This finding agree with the evidence of massive tidal tails 
around these two clusters (see above) indicating their strong 
interaction with the Galactic potential. For the other clus- 
ters of the sample, internal processes can be considered the 
main mechanisms of dynamical evolution. 

We revealed a slight correlation between 7 and the clus- 
ter mass and half-mass density which can be interpreted 
as a consequence of the fact that the most massive and 
dense tidally unaffecte d clusters are curre ntly in an expan- 
sion dominated phase (|Gieles et al.ll20l"ll ). The connection 
between slope 7 and half-mass density is less significant in 
the larger sample by MvM05 which is however limited to the 
innermost fit of the clusters. It is therefore not clear if such 
a relation is spurious (due to the small number of clusters of 
our sample) or real (and masked by the presence of outliers 
in the MvM05 sample). 
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